BOREL MEASURES WITH A DENSITY ON A COMPACT 
SEMI-ALGEBRAIC SET 
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JEAN B. LASSERRE 



Abstract. Let K C M" be a compact basic semi-algebraic set. We provide a 
necessary and sufficient condition (with no d priori bounding parameter) for 
a real sequence y = (ya), a £ N", to have a finite representing Borel measure 
absolutely continuous w.r.t. the Lebesgue measure on K, and with a density 
in n2^jLp(K). With an additional condition involving a bounding parameter, 
the condition is necessary and sufficient for existence of a density in Loo(K). 
Moreover, nonexistence of such a density can be detected by solving finitely 
many of a hierarchy of semidefinite programs. In particular, if the semidefinite 
program at step d of the hierarchy has no solution then the sequence cannot 
have a representing measure on K with a density in Lp(K) for any p > 2d. 



1. Introduction 

The famous Markov moment problem, (also called the L-problem of moments) is 
concerned with characterizing real sequences (s„), n S N, which are moment of a 
Borel probability measure ^ on [0, 1] with a bounded density with respect to (w.r.t.) 
the Lebesgue measure. It was posed by Markov and later solved by Hausdorff with 
the following necessary and sufficient condition: 

(1.1) l = so and < s^j < c/{n + 1), Vti, j e N, 

for some c > 



0, where s„j :— (—1)" ^(")A" ^Sj, and A is the forward operator 
in+i — Sn- Similarly, with p > 1, if in (jl.ll) one replaces the condition 



3nj 



< c/{n + 1) for all n,j e 



with the condition 



(1.2) 



\ 3=0 



1/p 



l)Sr. 



< c, 



then one obtains a characterization of real sequences having a representing Borel 
measure with a density in Lp([0, 1]) with p-norm bounded by c. For an illumi- 
nating discussion with historical remarks the reader is referred to Diaconis and 
Freedman [3] where the authors also make a connection with De Finetti's theorem 
on exchangeable 0-1 valued random variables. In addition, Putinar [9|[10] has pro- 
vided a characterization of extremal solutions of the two-dimensional L-problem of 
moments. 

An alternative if and only if characterization is via Linear Matrix Inequalities 
(LMIs) involving the moment and localizing matrices associated with the real se- 
quence and the Lebesgue measure A on [0,1], respectively. It states that the moment 
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and localizing matrices of the sequence must be dominated by that of A (scaled with 
some factor c > 0). This latter characterization carries over arbitrary compact ba- 
sic semi-algebraic sets of R" and a similar characterization also holds for measures 
on non-compact sets, satisfying an additional Carleman type condition; see e.g. 
Lasserre [3 p. 67-69]. 

A common feature of the above characterizations is the a priori unknown scalar 
c > which is precisely the required bound on the density. So for instance, if 
the conditions (|l.ip fail one does not know whether it is because the sequence has 
no representing measure fj, with a density in Loo([0, 1]) or because c is not large 
enough. So it is more appropriate to state that (jl.ll) are necessary and sufficient 
conditions for /i to have a density / e Loo{[Q, 1]) with ||/||oo < c. And similarly for 
(fTll for / e Lp([0, 1]) with II/IIp < c. 



Contribution: Consider a compact basic semi-algebraic set K C M" of the form 

(1.3) K:={xeM" : 5,(x) > 0, j = 1, . . . ,m}, 

for some polynomials gj G K.[x], j ^ 1, . . . ,m. For every, p S N, denote by Lp(K) 
the Lebesgue space of functions such that Jj^ |/|PA(dx) < cx). We then provide a 
set of conditions (S) with no a priori bound c, and such that: 

- A real sequence y = (j/a), a S N", has a finite representing Borel measure 
with a density in n2°^]^Lp(K), if and only if (S) is satisfied. In particular, if (S) 
is violated we obtain a condition (with no a priori bounding parameter c) for non 
existence of a density in npLp(K) (hence no density in _Loo(K) either). 

- A real sequence y = {ya), a G N", has a finite representing Borel measure with 
a density in Loo(K) if and only if (S) and an additional condition (involving an d 
priori bound c > 0), are satisfied. 

In addition, the conditions (S) consist of a hierarchy of Linear Matrix Inequalities 
(LMIs) and so can be tested numerically via available semidefinite programming 
softwares. In particular, if a finite Borel measure /i does not have a density in 
npLp(K) (and so no density in Loo(K) either), it can be detected by soh/ing finitely 
many semidefinite programs in the hierarchy until one has no feasible solution. That 
is, it can be detected from finitely many of its moments. This is illustrated on a 
simple example. 

Conversely, if the semidefinite program at step d of the hierarchy has no solution, 
then one may conclude that the real sequence y cannot have a representing measure 
on K with a density in Lp(K), for any p > 2d. 

So a distinguishing feature of our result is the absence of an d priori bound c in 
the set of condition (S) to test whether y has a density in n2^]^Lp(K). Crucial for 
our result is a representation of polynomials that are positive on K x R, by Powers 
[8]; see also Marshah [6l[7]. 

2. Main result 

2.1. Notation, definitions and preliminary results. LetM[x,i] (resp. M[x, t]d) 
denote the ring of real polynomials in the variables x = {xi, . . . ,Xn,t) (resp. poly- 
nomials of degree at most d), whereas 5][x, t] (resp. I][x, t]d) denotes its subset of 
sums of squares (s.o.s.) polynomials (resp. of s.o.s. of degree at most 2d). For 
every a S N" the notation x" stands for the monomial x^^ ■ ■ ■ x^" and for every 



NONNEGATIVITY 3 

deN, let N2+^ := {/3 G W'+^ ; J^j Pi < d} whose cardinal is s{d) = ("+;^+^). A 
polynomial / G M[x, t] is written 

(X,t)^/(x,t) = J2 fo^k^^t^ 

(Q,fc)eN"xN 

and / can be identified with its vector of coefficients f = {fak) in the canonical 
basis (x",*''), (a, /c) G N" x N, of IR[x, t]. But we can also write / as 

(2.1) (x,i)^/(x,t) = ^/fe(x)i^ 

fceN 
for finitely many polynomials f^ G M[x]. 

A real sequence z = (zak), {ct, k) G N" x N, has a representing measure if there 
exists some finite Borel measure i^ on M" x M. such that 

Zafc = / x"t''di/(x,t), V(a,fc)GN"xN. 

Given a real sequence z = (zq.^) define the linear functional Ly : ]R[x, f] — > M by: 

/(=^/afcx"t'=) ^iz(/) = 5I/"'«2/afe, /elR[x,i]. 

Moment matrix. The moment matrix associated with a sequence z = {z^k), 
{a, k) G N" x N, is the real symmetric matrix Mc;(z) with rows and columns indexed 
by N^"*""^, and whose entry (a, P) is just Za+p, for every a, /3 G N^"*"^- Alternatively, 
let Vd((x,i)) G »"('') be the vector ((x,i)"), a G N2+\ and define the matrices 

(B„) C S''^''^ by 



(2.2) Vd((x,i))v<i((x,t))^ = Yl B„(x,t)", V(x,t)G 



n+1 



Then Md(z) = I]QgN"+i ^cBc. 

If z has a representing measure i/ then Md(z) ^ because 



{i,Md{z)i) = f fdiy>0, Vf G 



K'^(rf). 



Localizing matrix. With z as above and 5 G IR.[x, t] (with (7(x, t) — ^ 37 (x, t)'''), 
the localizing matrix associated with z and g is the real symmetric matrix M(j(gz) 
with rows and columns indexed by N^"*" , and whose entry (a, /3) is just ^ g^Za+fs+^y , 
for every a,/3 G N;J+^ Ahernatively, let C G S''^'^^ be defined by: 

(2.3) 5(x,i)vd(x,t)vd(x,t)^ = ^ C„(x,i)", V(x,i)GM"+i. 

Then Md{gz) = J2aPN'^+^ ^a Cq. 

^ 2d+degg 

If z has a representing measure z^ whose support is contained in the set {(x, t) : 
5(x, t) > 0} then Md(5z) >z because 

{i,Md{9y)i) = I fgd,y>0, Vf G M^^'^^. 

With K as in (|1.3p . and for every j = 0, 1, . . . , m, let Vj := \{deggj)/2'] . 
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Definition 2.1. With K as m (El) let P{g) C M[x,i] be the 



convex cone: 



(2.4) Pig) = I Y. ^/3(x,t)5i(x)^^---5™(x)^" : V^eS[x,t] 

[ /36{0,1}" 

The convex cone P{g) is called a preordering associated with the (?j's. 

Proposition 2.2. Le^ K be as in \1.!^) . A polynomial f G M[x, t] is nonnegative 
on K X M only if f can be written as 

2d 

(2.5) (x,t) ^^ /(x,t) = ^/,(x)i^ 

fc=0 

/or some rf G N anrf where f2d > on K. 

Proof. Suppose that the highest degree in t is 2d+l for some d G N. Then /2d+i 7^ 
and so by fixing an arbitrary Xq G K, the univariate polynomial t 1-^ /(xo,i) can 
be made negative, in contradiction with / > on K x R. Hence the highest degree 
in t is even, say 2d. But then of course, for obvious reasons f2d > on K. D 

We have the following important preliminary result. 

Theorem 2.3 (011]). Let K as in lll.3\} be compact and let f G R[x, i] be of the 
form f{x,t) — "^^=0 ■fk{^)t'^ for some polynomials (f^) C IR[x], and with f2d > 
on K. Then f G Pig) if f >0 onKxR. 

And so we can derive a version of the K x M-moment problem where for each 
/? G N™, the notation g^ stands for the polynomial g^ ■ ■ ■ g!^". 

Corollary 2.4. Let K as in il.3\} be compact. A real sequence z = izak), (a, fc) G 
N" X N. has a representing measure on K x M. if and only if 

(2.6) Md(z)^0; Mdig^z)hO, /3 G {0, 1}™, 

for every d G N. 

Proof. The only if part is straightforward from the definition of the moment and 
localizing matrix M(j(z) and NLdig^ z), respectively. 

The if part. Suppose that (|2.6|) holds true, and let / G M[x, t] be nonnegative 
on the closed set K x M. Hence by Proposition 12. 2) / has the decomposition ()2.5p 
for some integer d 7^ 0. For every e > 0, the polynomial (x, t) 1-^ /^(x, i) := 
/(x,i) + e(l + t^'') has the decomposition 

2d 



/,(x,i) = ^/,fc(x)t^ 



fe=0 

with /eo = /o + £ ^nd /e2d(x) = /2d(x) + £. Therefore, /e is strictly positive on 
K X K, and /,2d > on K. By Theorem[231 /, G Q(g), i.e., 

/e(x,i) = ^ VMx,t)ff(x)^ 

/3e{o,i}™ 
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for some SOS polynomials {^pp) C S[x,i]. Next, let z satisfy ()2.6|) . Then 

L,(/)+eLy(l+i2'i) = i,(/^) 

= ^ iz(#/)>0 
/3e{o,i}- 

where the last inequality follows from 

Md{g^ z) ^ ^ L,{h^ /) > 0, V/i £ R[x, t]d, 

for every /? G {0, 1}™. But since Lz(l + t'^'^) > and e > was arbitrary, one 
may conclude that L^{f) > for every / e M[x, t] which is nonnegative on K x R. 
Hence by the Riesz-Haviland theorem (see e.g. ^ Theorem 3.1, p. 53]), z has a 
representing measure on K x M. D 

2.2. Main result. Let Loo(K) be the Lebesgue space of integrable functions on 
K (with respect to the Lebesgue measure A on K, scaled to a probability measure) 
and essentially bounded on K. And with 1 < p < oo, let Lp(K) be the Lebesgue 
space of integrable functions / on K such that /^ |/|^A(dx) < oo. A Borel measure 
fj, absolutely continuous w.r.t. A is denoted /i <C A. 

Theorem 2.5. Let K C M" &e as m f73)j and let y = (y^), a e N", be a 
real sequence with yo ~ 1. Then the following two propositions (i) and (ii) are 
equivalent: 

(i) y has a representing Borel probability measure /i <S A on K, with a density 
m n^iLp(K). 

(ii) Mrf(y) > and Mrf(/ y) ^ for all (5 e {0, 1}" and all deE. In addition, 
there exists a sequence z = (zak), (o^, k) G N" x N, such that i2. 6\) holds, and 



Zao = / x" A(dx) ; 



(2.7) z„o = / x" A(dx) ; z„i = j/,, Va G N". 



Moreover, if in |^. 7| ) one includes the additional condition supj, zok < cxd, then (ii) 
is necessary and sufficient for y to have a representing Borel probability measure 
/I <C A on K, with a density in Loo(K). 

Proof. The (i) => (ii) implication. As y has a representing Borel probability mea- 
sure /i on K with a density / G ip(K) for every p = 1, 2, . . ., one may write 



fiiA) = / /(x)A(dx), VA G 

JA 

Define the stochastic kernel (p{B\x.), B G B{R), x G K, where for almost all x G K, 
(/?(• I x) is the Dirac measure at the point /(x). Next, let v be the finite Borel 
measure on K x R defined by 

(2.8) u{A xB) -.^ [ (p(B|x) A(dx), VA G i3(R"), B G B(R). 

J A 

Let z = (yak), (a, fc) G N" x N, be the sequence of moments of i'. 

Zak = [ x°'t''diy{x,t) = /" x" ( f t''tp{dt\x)) \{dx), 
Jr Jk \Jr J 

(2.9) = / x" /(x)*^ A(dx) (weU defined as / G Lp(K) for ah p). 

JK 
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In particular, for every a G N", 

Zao = / x" A(dx); Zai = / x"/(x)A(dx) = / x" d^ == y^. 
Jk Jk Jk 

Moreover, as ly is supported on K x R then Mc;(y) ^ and 'M.dig'^ y) ^ for all d 
and all /3 e {0, 1}™. Hence dSHl-dHZl) hold. 

The (ii) ^ (i) imphcation. Let z = (z„fc) be such that ([^ - (PT)) hold. By 
Corollary (|2.4p . z has a representing Borel probability measure j/ on K x R. One 
may disintegrate v in the form 

v{AxB)=( tp{B\x)'tP{dx), B eBiR), AeB{R"), 
JAnK 

for some stochastic kernel ip(- \ x), and where ip is the marginal (probability mea- 
sure) of f on K. From (|2.7p we deduce that 



/ x"V(rfx) = z„o = / x"A(dx), 
JK Jk 



Va e N", 
which, as K is compact, implies that t/j ~ A. In addition, still from ()2.7p . 

z„i = / x"idz/(x,i) = /" x" ( /"t(p(di|x)j A(rfx) VaeN" 

v ' 

/(x) 

(2.10) = / x" /(x) A(dx) VaeN", 

de(x) 

where / : K — > R is the measurable function x H> Jg^tip{dt\x), and 6* is the signed 
Borel measure 9{B) := J^^^ /(x) A(dx), for all B e B(R"). 

But as K is compact, by Schmiidgen's Positivstellensatz |TT1, the conditions 

Md(y)^0, Md(/y)^0, /3 e {0, 1}™, Vd G N, 

imply that y has a finite representing Borel measure /i on K. And so as Zai = ya 
for all a G N", and measures on compact sets are moment determinate, one may 
conclude that dfi = d9 ~ f d\, that is, /i ^ A on K. Next, observe that for every 
pGN, 

zop = / tPdiy{y:,t) ^ f ( f t^ ip{dt\xy\ A(dx) Va G N" 



/(x)p 
/(x)PA(dx) VaGN", 

'K 

and so / G ip(K) for all p > 1. 

Finally consider (|2.7p with the additional condition sup„ zop < cx). Then in the 
above proof of (i) =^ (ii) and since now y has a finite representing Borel measure 
with a density / G Loo(K), one has limp_5.oo ||/||p — ||/||oo because K is compact; 
see e.g. Ash [51 problem 9, p. 91]. And therefore since zop = /k /(^)^^('^^)' "^^ 
obtain supp zop < oo. 

Similarly, in the above proof of (ii) ^ (i), sup^ ZQp < oo implies sup^ Jj^ /(x)PA((ix) 
supp II /lip < oo. But this implies that / G ioo(K) since K is compact. D 
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Computational procedure. Let 7 = (7a), ct G N", the moment of the Lebesgue 
measure on K, scaled to make it a probabihty measure. In fact, the (scaled) 
Lebesgue measure on any box that contains K is fine. 

Let y = (j/q), a G N", be a real given sequence, and with K as in (|1.3p let 
Vj := \{deggj)/2], j ~ 1, . . . ,m. To check the conditions in Theorem 12.5^ 11). one 
solves the hierarchy of optimization problems, parametrized by d e N. 

Pd = min trace(M(;(z)) 

s.t. Md(z) h 
(2-11) Md_„^(/z) h 0, /3e{0,l}™ 

Zao = 7q: a e N^^ 
2ai = 2/a, (a,l) G^M^- 

Each problem (|2.1ip is a semidefinite progranu. Moreover, if (|2.1ip has a feasible 
solution then it has an optimal solution. This is because as one minimizes the trace 
of Md(z), the feasible set is bounded and closed, hence compact. 

In p. lip one may also include the additional constraints zq^ < c, k < 2d, for 
some fixed c > 0. Then by Theorem 12.51 y has a representing Borel probability 
measure on K with a density in Loo(K) bounded by c, if and only ii pd < 00 for all 
d. 

Each semidefinite program of the hierarchy (|2.1ip . cf G N, has a dual which is 
also a semidefinite program and which reads: 
(2.12) 

p*,— iasiyi I p(x)A(dx) + iy(g) 

m 

S.t. Y. (x"t'=)'-(p(x)+tg(x)) = ao(x,i) + ^a,(x,i)g,(x) 

(Q,fc)GN2 + ^ i=l 

degp< 2d; degq < 2rf- 1; Gj G S[x, i]f_„^, j = 0,...,m, 

where uq = 0. In particular, if y is the sequence of a Borel measure on K then in 

(12.121) one may replace Ly{q) with / q{x.)dfi{-x). 

Jk 

2.3. On membership in Lp(K). An interesting feature of the hierarchy of semi- 
definite programs (|2.11l) , d G N, is that it can be used to detect if a given sequence 
y = iVa), 01 G N", cannot have a representing Borel measure on K with a density 
in £p(K), p > 1. 

Corollary 2.6. Let K C K" he as m IfT^) and let y = (y„), a G N", be a 
real sequence with yo = 1. If the semidefinite program 112.11]) with d G N, has no 
solution then y cannot have a representing finite Borel measure on K with a density 
in LpCK.), for any p > 2d. 

Proof. Suppose that y has a representing measure on K with a density / G L2d(K), 
and hence in ifc(K) for all k < 2d. Proceeding as in the proof of Theorem 12. 5[ let 
f be the Borel measure on K x M defined in (12.81). Then from (12.91) one obtains 



zck = / x" /(x)*'- A(dx), (a, k) G N^+\ 



A semidefinite program is a convex optimization problem that can be solved efficiently, i.e., 
up to arbitrary fixed precision it can be solved in time polynomial in the input size of the problem; 
see e.g. [l]. 
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which is well-defined since K is compact (so that x" is bounded) and k < 2d. And 
so the sequence z — (zafc), («, k) S Njj"^, is a feasible solution of (|2.11l) with d. D 

Notice that again, the detection of absence of a density in ip(K) is possible with 
no a priori bounding parameter c. But of course, the condition is only sufficient. 

Example 1. Let K := [0,1] and s € [0,1]. Let A be the Lebesgue measure on 
[0, 1] and let 6s be the Dirac measure at s. One wants to detect that the Borel 
probability measure fj.a '■= aA + (1 — a)6s, with a G (0, 1) has no density in Loo(K). 
Then (1213 reads 



Zko = , -. , fc = 0, 1, . . . ; Zki = -r—— + (1 - a)s'', fc = 0, 1, . . . 

The set K is defined by {x : g{x) > 0} with x H> g{x) := x{l — x). We have tested 
the conditions Mc;(z) ^ and Md(gz) ^ along with (|2.7p where k < 2d (for zpfc) 
and k <2d— 1 (for z^i ) . 

We have considered a Dirac at the points s — fc/10, k = 1, . . . , 10, and with 
weights a = l — fc/10,fc = l,...,10. To solve (I2.11[) we have used the GloptiPoly 
software of Henrion et al. [3] dedicated to solving there generalized problem of 
moments. Results are displayed in Table [T] which should be read as follows: 

• A column is parametrized by the number of moments involved in the conditions 
(|2.7p . For instance. Column "10" refers to (|2.7I) with d — 10/2, that is, the moment 
matrix Md(z) involves moments Zij with i + j < 10, i.e., moments up to order 10. 

• Each row is indexed by the location of the Dirac Sg, s G [0, 1] (with /i^ — 
oA + (1 - a)Ss). The statement "1 - a > 0.5" in row "s = 0.3" and column "10" 
means that (12.71) is violated whenever 1 — a > 0.5, i.e., when the weight associated 
to the Dirac Ss is larger than 0.5. 

One may see that no matter where the point s is located in the interval [0, 1] , if its 
weight 1 — a is above 0.5 then detection of impossibility of a density in Loo([0, Ij) 
occurs with moments up to order 10. If its weight 1 — a is only above 0.1 then 
detection of impossibility occurs with moments up to order 12. So even with a 
small weight on the Dirac Sg , detection of impossibility does not require moments 
of order larger than 12. 

Example 2. Still with K = [0, 1], consider now the case where fia = aX+{l—a){Ss+ 
(5s+o.i)/2, that is, fia is a (a, 1 — a) convex combination of the uniform probability 
distribution on [0, 1] with two Dirac measures at the points s and s + 0.1 of [0, 1], 
with equal weights. The results displayed in Table [2] are qualitatively very similar 
to the results in Table [1] for the case of one Dirac. 
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s \ moments 


8 


10 


12 


14 


0.0 


1 - a > 0.3 




-a > 0.1 


l-a> 0.1 




-a > 0.1 


0.1 


l-a> 1 


1 - 


-a > 0.3 


1 -a > 0.2 




-a > 0.2 


0.2 


l-a> 1 


1 - 


-a > 0.3 


l-a> 0.1 




-a >0.1 


0.3 


l-a> 1 


1 - 


-a > 0.5 


l-a> 0.2 




-a >0.1 


0.4 


l-a> 1 


1 - 


-a > 0.5 


l-o> 0.2 




-a >0.1 


0.5 


l-a> 1 


1 - 


-a > 0.5 


l-a> 0.2 




-a >0.1 


0.6 


l-a> 1 


1 - 


-a > 0.5 


l-a> 0.2 




-a >0.1 


0.7 


l-a> 1 


1 - 


-a > 0.5 


l-o> 0.1 




-a >0.1 


0.8 


l-a> 1 


1 - 


-a>0.5 


l-a> 0.1 




-a > 0.1 


0.9 


l-a> 1 


1 - 


-a>0.5 


1-0 7^0.4,0.5 




-a > 0.1 


1.0 


1 - a > 0.4 


1 - 


-a> 0.1 


l-a> 0.1 




-a >0.1 



Table 1 . Moments required for detection of failure; one Dirac 



(s, s + 0.1) \ moments 


10 


12 


(0.1,0.2) 


1 - a > 0.4 


l-a> 0.2 


(0.2,0.3) 


1 - a > 0.6 


l-a> 0.2 


(0.3,0.4) 


1 - o > 0.5 


l-a> 0.2 


(0.4,0.5) 


1 - a > 0.6 


l-a> 0.1 


(0.5,0.6) 


1 - a > 0.6 


1 - a > 0.2 


(0.6,0.7) 


l-a>0.7 
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